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Double parton scattering is sensitive to correlations between the two partons in the hadron, 
including correlations in flavor, spin, color, momentum fractions and transverse separation. We 
obtain a first estimate of the size of these correlations by calculating the corresponding double 
parton distribution functions in a bag model of the proton. We find significant correlations between 
momentum fractions, spin and flavor, but negligible correlations with transverse separation. The 
model estimates of the relative importance of these correlations will help experimental studies 
disentangle them. 



I. INTRODUCTION 

High energy scattering processes such as Drell-Yan pro- 
duction, pp — > £ + £~, are described by the scattering of 
two incoming partons, and the cross section is given by 
the convolution of a partonic scattering cross section a 
and parton distribution functions (PDFs). Sometimes 
two hard partonic collisions take place within a single 
hadronic collision, a process which is known as double 
parton scattering (DPS). DPS is higher twist, i.e. it is 
suppressed by a power of Aq CD /Q 2 , where Q is the par- 
tonic center-of-mass energy of the collision. Despite this 
power suppression, the DPS scattering rate is still large 
enough that it has become a background for new physics 
searches at the LHC. For example, DPS contributes to 
same-sign WW and same-sign dilepton production [I]- 
4 , and is a background for Higgs studies in the channel 
pp -> WH — s- ivbb [SHS] . DPS has been observed at the 
LHC; a preliminary study using 33 pb -1 of data found 
that 16% of the W + 2 jet events were due to DPS 0. 

In the original work on DPS, the cross section was 
written as [10] 

do- = i ^2 Jd 2 z ± F l: j(x 1 ,X2,z ± ,fi)F k i(x 3 ,X4 : ,z_ L ,n) 

x a lk (xixz\fs, n)a j i(x 2 x i ^,n) . (1) 

The double parton distribution function (dPDF) 
F i j(xi,x 2 ,z±) describes the probability of finding two 
partons with flavors i, j = g,u,u,d, . . . , longitudinal mo- 
mentum fractions x%, x 2 and transverse separation zj_ in- 
side the hadron. The partonic cross sections a describe 
the short distance processes, and S is a symmetry fac- 
tor that arises for identical particles in the final state. 
Eq. (JlJ ignores additional contributions that are sensi- 
tive to diparton correlations in flavor, spin and color, as 
well as parton-exchange interference contributions 
03]- These correlations are present in QCD, and one of 
our goals is to estimate the size of these effects. 

It is commonly assumed in DPS studies that the depen- 
dence on the transverse separation is uncorrelated with 



the momentum fractions or parton flavors, 

Fij(xi,x 2 ,z_L,IJ>) = Fij(xi,X2,n)G(z ± ,n) . (2) 

In addition, a factorized ansatz is often made 

Fij(xi,X2,fjt) = fi(xi,n)fj(x 2 ,n) 6(l-x 1 -x 2 ){l-x 1 -x 2 ) 

(3) 

where / denotes the usual (single) PDF. The factor 
9(1 — x\ — x 2 ){l — x\ — x 2 ) n smoothly imposes the kine- 
matic constraint x\ + x 2 < 1, and different values of the 
parameter n > have been considered. 

The dPDF is a nonperturbative function, but once it is 
known at a certain scale /i, its renormalization group evo- 
lution can be used to evaluate it at a different scale. The 
evolution of F^jixx, x 2 , A*) was determined a long time 
ago [Til 023- It has recently been extended to include 
the zj_ dependence and describe correlation and interfer- 
ence dPDFs [HHl"4l IT7] . The color-correlated and inter- 
ference dPDFs are all Sudakov suppressed at high ener- 
gies Q31 HB] and will therefore not be considered. 

Eventually the dPDFs will be determined by fitting to 
experimental data, just as for the usual PDFs. Ref. [19] 
goes a step in this direction, showing how angular cor- 
relations in double Drell-Yan production may be used to 
study spin correlations in dPDFs. In this paper, we de- 
termine the dPDFs at a low scale /i ~ Aqcd using a bag 
model for the proton [5D]. This model calculation pro- 
vides an estimate of the importance of various diparton 
correlations, which can be used to guide the experimen- 
tal analysis. It also provides an estimate of dPDF distri- 
butions in the absence of more accurate determinations 
directly from experiment. 

We follow some of the existing structure function cal- 
culations in the bag model [2"TM23| . There are obvious 
limitations to this approach, just as for bag model cal- 
culations of the usual PDFs. First of all, the bag model 
only describes valence quarks. Bag model calculations 
are only meaningful when the fields in the dPDF act 
inside the bag, which restricts the momentum fractions 
x > 1/(2MR) ^0.1, where M is the proton mass and R 
is the bag radius. Finally, the bag was treated as rigid 
in the early literature, Ref. [21] . A consequence is that 
momentum is not conserved, and parton distributions do 
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not vanish outside the physical region (x > 1). Alterna- 
tive treatments of the bag were proposed to alleviate this 
problem [22 24] . We emphasize that we are not attempt- 
ing to use the most sophisticated bag model description 
of the proton. Rather, we simply want to provide a first 
estimate of the size of the various correlation effects. Bag 
model PDFs are usually chosen as the initial value of 
PDFs at a low scale \i ~ Aqcd, which are then evolved 
to higher scales using their QCD evolution. Since in the 
bag model the valence quarks carry all the momentum, 
this initial scale fi needs to be taken quite low [2"2~| . 

We also investigate Eqs. ^ and ([3| in this paper, us- 
ing our bag model results. We find that Eq. ^ holds 
reasonably well, but Eq. (|3| is badly violated. Problems 
with Eq. ([3]) have already been pointed out in Ref. [25] . 
using sum-rules and the evolution of the dPDF. In the 
simplest bag models of the type we consider, the color- 
correlated dPDFs F T are given by —2/3 times the color- 
direct dPDFs F 1 , since diquarks are in a 3 representation 
of color. 



II. CALCULATION 

We briefly summarize the ingredients of the bag 
model [20] that are needed to calculate the dPDFs. The 
bag model wave functions are the solutions of the mass- 
less Dirac equation in a spherical cavity of radius R. We 
only need the ground state, which is given by 



* m (r,t) = N e 



iflt/R 



for a bag centered at the origin. Here f2 = ER, 



N 2 = 



1 



n 4 



r 3 n 2 - sin 2 n ' 



47T V° m ^ 



(4) 



(5) 



ji are spherical Bessel functions, and rn =t,|- The con- 
dition that the color current does not flow through the 
boundary r^j^T^ \\ t \=r = leads to 



j (n)=h(Q) 



Q fa 2.043 , 



(6) 



and we will take R = 1.6 fm in our numerical analysis. 

The quark field is expanded in terms of bag wave func- 
tions, quark creation and annihilation operators £lj(a), 
aj(a) and antiquark creation and annihilation operators 
fci(a), b\(a). These operators create or annihilate quarks 
and antiquarks in a bag centered at r = a [sec Eq. (14)]. 
We only need the operators for the ground state wave 
functions in Eq. Q with m =f, i- 

The spin-up proton wave function is given in terms of 
the standard quark model wave functions as 



V6 



uud) (2im> - int> - ntt» 



(I) 



As usual, the color indices are suppressed, and the wave 
function has to be symmetrized over permutations. Ig- 
noring color, one can also write the wave function in 
terms of bosonic [2S] creation operators, 



|Pt,r = a) = 



1 



4t( a )4t( a )4i( a ) 



V3 



|0,r 



(8) 



Here \P f, r = a) and |0, r = a) are the proton and empty 
bag state, respectively, both at position a. The a£ m (a) 
denotes the creation operator for a quark of flavor q with 
spin m in a bag at position a. 

An important difference between various calculations 
in the literature is the treatment of the overlap between 
empty bags at different positions, 

(0,r = a|0,r = b) =5 3 (a-b) in Ref. HJ , 

(0,r = a|0,r = b) = 1 in Refe. [21 [23] . (9) 

These opposite limits treat the bags as either completely 
rigid or fully flexible, and the latter will be our default. 
We will return to the rigid bag in Sec. |IID| To account 
for the displacement between bags, we follow Ref. [22] in 
taking 

{oi(a),ot(b)} = % J d 3 a;*t(x - b)^(x - a) . (10) 

For the rigid bag these are replaced by the familiar anti- 
commutation relations 

{a i ,a]} = S ij , (11) 

where we only need the relation when a and cl* are at the 
same bag position, because of Eq. ([9| . 

The proton state with momentum p is constructed us- 
ing the Peierls-Yoccoz (PY) projection [2"T] . 

\P,p) = ~ /d 3 a e ia 'P|P,r = a), (12) 

03 (P) J 

where <ps(p) fixes the (non-relativistic) normalization of 
the state. The functions 0„(p) are given by 



dae" ip - a /dx* t (x-a)*(x) , (13) 



\Mp)\' 



which wc will need for n = 1, 2, 3. 

The final ingredient is the expression for quark fields 
acting in the bag. The field for a u-quark relative to a 
bag at a is given by [22] 



m(x, t) = ^ a «m(a)* TO (x-a)e 1 



Qt/R 



(14) 



Here ". . . " denotes contributions from other bag states 
that will not be needecQ The expression for d-quarks is 
similar. 



1 We will not consider the so-called z-graph or four-quark inter- 
mediate state contribution |21l 1231 . where the field creates an 
antiquark. This only contributes at small x and is thus outside 
the range of validity of the calculation. 
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function sets 



FIG. 1. The proton PDFs u (solid red), Au (dashed blue) 
and Su (dotted green). 



A. Single PDF 

We first summarize the well-known calculation of the 
(single) PDF in the bag model. The light-cone vectors 
are 
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= (1,0,0,1), = (1,0,0,-1) . 



(15) 



and we assume the light-cone gauge n ■ A = 0. In the 
proton rest frame, where p M = (M, 0)J^] 



q(x) = 2M 



dz + 



-ixMz + /2 



47T 

x (p,p = o|g( 



(16) 



Pp = 



(P,r = 0|at roi (0)a 9TO2 (0)|P,r = 0) 

mi,tH2=t4 

2M d MW *„ lW f*„ W l|« 



(2„)» 



•\n/R-xM\ 

The factor of 2M is due to the nonrelativistic normal- 



v- 



R 



xM . 



(17) 



In using Eq. (14) we assumed that the field q acts at 
the position of the bag of the left state and q at the 
position of the bag of the right state [25] . The third line 



of Eq. ( 16 ) contains all the dependence on the spin-flavor 
wave function of Eq. ^ , which is connected with the spin 
of the bag wave functions through the sum on m\^. For 
the unpolarized PDF only m\ — 1112 contributes, and the 
matrix element simply counts the number of quarks of a 
given flavor q in the proton, 



= E <^ = 0| 



qm 



(0) a<?m (0)|P,r = 0). (18) 



An immediate implication of Eq. ( 16 ) is that the peak of 
q is located at x = fi /(MR), independent of the quark 
flavor. This disagreement with experimental measure- 
ments of u and d may be alleviated by refining the bag 
model, see e.g. Ref. [35]. We will restrict ourselves to 
the simplest bag models in this paper, so its limitations 
should be be kept in mind while using the results. 

The extension of Eq. (16) to longitudinal and trans- 



ization of states, and 02 is given by Eq. (13). The delta 



versely polarized quark distributions is given by replac- 
ing I by ^75 for Aq and §7^75 for 5q. Aq and 5q 
only contribute in processes involving longitudinally and 
transversely polarized protons, respectively. The matrix 
elements required are evaluated in Sec. |II C| To aid the 
evaluation of the remaining integral in Eq. (16), conve- 
nient expressions for the functions <f>i and the bag wave 
function in momentum space are given in Appendix [A] 
The resulting distributions for the various proton PDFs 
are compared in Fig. [T] 



B. Double PDF 

We calculate the double PDF using the definitions in 
Ref. [14]. We will not consider color correlated or inter- 
ference double PDF, since these are Sudakov suppressed. 
Starting with the spin-averaged dPDF Fg ig2 (xi,X2,kj_) 



(xi,x 2 ,kj_) = -8ttM 2 /dVe" iZi ' kj 



dz 1 dz 2 dz 3 c -j Xl Mzt /2 e -ix 2 Mz£ /2 ^\x\Mz^ /2 

Air 47r 47r 



(19) 



x P,p = 



<h 



zX — + zj_ 



P,p = 



E ( P ' r = °\ a l** (0)<4m 2 Wa q2m4 (0)a qima (0) | P, r = 0) 



mi ,ni2 ,7713,3714 



2 We also use the notation q for the PDF f q , and qq, qAq, . . . for 
the dPDFs F qq , F qAq , .... 
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x (27r) 2 5 2 (k 1± - k 3± - k x ) * mi (k x )f * m3 (k 3 ) * m3 (k 2 )f * m4 (k x + k 2 - k 3 ) l01 j kl + k 2 : 

2 2 |</>3(0)r 



r 



where </>i is given by Eq. ( 13 1. Results for the matrix ele- and we find the weighting: 



ments on the third line of Eq 



19 1 are given in Sec. II C 



The expressions in Appcndix[A aid an efficient evaluation 
of the remaining integrals. The numerical integrals were 
performed using the CUBA integration package |29j . 



C. Spin Correlations 

The computation of spin-correlated double PDFs is al- 
most identical to Eq. (fl9| . The spinors f ® § in Eq. (fl9|) 



are replaced by [T2"Hl"4] 
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We have multiplied the dPDFs by appropriate factors of 
M|zx| to remove explicit dependence on \z± \ and M from 
the spin structure. These will be the quantities that we 
plot. As in Eq. (19), we again switch to Fourier space. 



Azimuthal symmetry implies that we may simply pick 
kj^ to be along the 1 direction, resulting in the 7-matrix 
structures: 



M\z ± \F qi5q2 
M\z± \ F AqiSq2 

FSq 1 Sq 2 

X ^SqiSq 2 



ft 1 

f7 X 75 ® 57*75 



I ®h 2 
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C 
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The spin structures q\8q 2 and A<ji<5g2 vanish in our cal- 
culation. This follows from the factorized form of the 
dPDF in Eq. (19 1, where only 0i(ki + k2) introduces 
correlations between ki and k2. Specifically, this allows 
one to perform different reflections on each momentum, 
k lx -> |kx| - k lx and k 2x -> — |kj_| - k 2x , under which 
the integrand is odd. Though this is due to the form of 
the bag model matrix elements, it suggests that the spin 
structures q±Sq 2 and Aqi$q2 are likely smaller than the 
others. 

We now evaluate the spin-flavor matrix elements that 
enter in the single and double PDFs. Since we suppressed 
the antisymmetric color wave function of the proton, 
the creation and annihilation operators essentially satisfy 
commutation relations. For the unpolarized and longitu- 
dinally polarized single PDF only 777,1 = 7712 contributes, 
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For Sq we need a transversely polarized proton 
1 
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(20) 



The non- vanishing matrix elements are 
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The double PDFs we consider are invariant under spin 
flip (they are only sensitive to diparton spin correlations), 
so we can simply use a spin-up proton. The nonvanishing 
matrix elements are 
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Note that due to these spin-flavor correlations, the dPDF 
for uu and ud do not simply differ by an overall factor, 
as is the case for the single PDF. 



D. Rigid Bag 

For a rigid bag, the overlap of empty bag states is 

(0,r = a|0,r = b) = <S 3 (a-b). (21) 
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FIG. 2. Plot of the PY factors which enter the calculation 
of the single PDF (dotted blue) and double PDF (solid red). 
They suppress the PDFs in the unphysical regions x > 1 and 
x < 0. 




FIG. 3. Plot of the bag model proton PDF u(x) with (solid 
red) and without (dotted blue) PY factors. 



The only change to the single PDF in Eq. ( 16 ) is that 
it removes the PY factor 1 2 (k) | 2 /|0 3 (0) | 2 . This factor 
suppresses the "leakage" of the PDF into the unphysical 
regions x < and x > 1, without affecting the integral 
over all x, see Sec. |IIE| The PY factor is plotted in Fig. [2] 
and the PDF with and without the PY factor is shown 
in Fig. [J 



Similarly, the rigid bag overlap in Eq. (21) removes 
the factor |0i(k)| 2 /|</> 3 (O)| 2 (also plotted in Fig. from 
the double PDF in Eq. Q. In this case the double PDF 
factors and there are no correlations between the momen- 
tum fractions x\ and x 2 , which is a clear shortcoming of 
treating the bag as rigid. At kj^ = the rigid bag dPDF 
takes a particularly simple form 



F qiq2 (x 1 ,x 2 ,'k_ L = 0) 



n qi n qi 



qi(x 1 )q 2 (x 2 ) , (22) 



where the coefficient c qiq2 is fixed by the spin-flavor wave 
function 



J <?1<?2 
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mi —77*3 ,7712—777-4 



(P,r = 0|aJ imi (0)4 2 m 2 (0) 



xa 92m4 (0)a gim3 (0)|P,r = 0). (23) 
From the tables in Sec. II C| we find that c uu = c u d = 2. 



E. Normalization 



The normalization of the single PDF and double PDF 
is given by integrating over all x, including unphysical 
regions. Both treatments of the bag in Eq. ^ will be 
considered. The single PDF in a rigid bag gives 
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Here we used that 
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i . 



(24) 



(25) 



This second equation and the symmetry between n and 
n implies that we could replace / 4 — >■ 1/2 in Eq. (24 1. 



The corresponding calculation with a flexible bag, i.e. 
including the PY factor, is 
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(26) 



and has the same normalization. However, the PY factor 
reduces the PDF at unphysical x. Specifically, 2% of 
the contribution to the integral in Eq. ( 26 1 is from the 



unphysical region, compared to 11% in Eq. (24) 



For the double PDF, the normalization for the rigid 
bag follows from Eqs. p2l and ((24 ) 
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(27) 
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FIG. 4. The double PDF uu{x\, X2, kx) as a function of x\ and |kj_| for fixed xi — 0.4. The right panel tests the ansatz in 
Eq. |2| that Xj and kx are uncorrelated. This holds reasonably well, since the different |kj_| curves are nearly identical. 





FIG. 5. The double PDF uu(xi, X2, kx) as a function of xi and X2 for fixed kx = 0. In the right panel we divide by u(x2) to 
test the often used assumption in Eq. (J3J that the Xi are uncorrelated. This clearly fails, since the ratio depends strongly on 
x 2 - 



where the coefficient 



-9l<?2 



is given in Eq. ( 23 1 . The cal 



culation including the PY factor is similar to Eq. ( 26 ) 
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(28) 



The small correction with respect to Eq. (27 1 arises be- 



cause we can no longer replace ?i^/4 — ¥ 1/2. Specifically, 
Eq. (|Al| implies 
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(fl 2 - sin 2 fl) ' 

Since the momenta k lz and k 2z become correlated 



through tfii (ki 
(k l2 )(k 2z ) = 0. 



k 2 ), this implies that (ki 2 k 2z ) i= 



III. PARTON CORRELATIONS 

We are now ready to investigate the size of the var- 
ious diparton correlation effects using the bag model 
dPDFs. We start by studying the dependence of the 
dPDF uu(xi, x 2 , kj_) on x\ and |kj_|, keeping x 2 = 0.4 
fixed for simplicity. As the left panel of Fig. [4] shows, 
the dPDF reduces significantly with increasing |kjj. In 
the right panel we test the ansatz in Eq. ([2| that the 
dependence on Xi and kx is uncorrelated, by dividing by 
uu(0A, 0.4, kj_). If the ansatz holds, the universal trans- 
verse function G(kjJ should drop out in this ratio, mak- 
ing the result independent of kj_ ■ As the plot shows, this 
seems to holds quite well. It only breaks down for the 
largest values of |kj_|, where the dPDF is orders of mag- 
nitude smaller than at |kjj =0. We also note that there 
is some leakage into the unphysical region x\ + x 2 > 1, 
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as was the case for the single PDF in Fig. [3j though this 
effect is reasonably small. 

Next we explore the x\,x 2 dependence of dPDF 
uu(x\, X2, kj_) for kj^ = 0, which is shown in Fig. [5} As 
x 2 is increased, the peak of the jq-distribution moves to 
smaller xi, responding to the reduced momentum avail- 
able. The peak height reduces as well, though not for 
small x 2 since the bag model only describes the valence 
quarks. To test the factorization ansatz in Eq. ^ for 
ri = 0, we divide by u{x 2 ) in the right panel. Since the 
resulting distributions clearly still depend on x%, corre- 
lations between X\ and x 2 are important. Inclusion of 
the factor of (1 — X\ — X2)™ does not alter this conclusion. 
The correlations can also be seen in the three-dimensional 
plot of Fig. (6) We remind the reader that this conclusion 
depends on the treatment of the bag, since X\ and x 2 
wou ld be uncorrelated if a rigid bag was assumed (see 

Sec.[nT3|). 

The relative size of the various spin structures in 



Sec. II C| are studied in Fig. [7] They are shown as a func- 
tion of X\ (top row) and |kj_| (bottom row), keeping all 
other variables fixed. All spin structures show a similar 
dependence on x\, though there is a hierarchy between 
their sizes. The shape of the |kj_| dependence differs be- 
tween the various spin structures. They all decrease for 
large |kj_|, but SqSq* also decreases as |kj_| goes to zero 
whereas the others are maximal at |kjj = 0. This can 
be understood because SqSq* explicitly requires a trans- 
verse scale (z± or kjj to be defined, so it must vanish 
at |kj_| =0. These same figures also illustrate the differ- 
ences between the uu (left column) and ud (right column) 
dPDF. Unlike the single PDF, where the difference be- 
tween u and d was simply an overall factor oin u /rid = 2, 
the double PDF has more flavor dependence. This arises 
through the spin dependence and the correlations in the 
spin-flavor wave function. As Fig.[7]shows, the difference 

UU(X] ,X2,k_ L = 0)/[u(X])u(X2)] 




FIG. 6. The correlation between the momentum fractions of 
two u quarks in the proton is shown by plotting the ratio 
of the double PDF uu(xi, xi, kj_ = 0) to the product of two 
single PDFs u(xi)u(x 2 ). 



between uu and ud is fairly small. However, the spin 
correlations are about twice as big for ud than for uu. 



IV. CONCLUSIONS 

We have computed the dPDFs using a bag model for 
the proton. The bag model results should be treated as 
the dPDFs at a low scale, which can then be evolved 
to higher energy using the known QCD evolution equa- 
tions [E1[T7]. We find substantial diparton correlations 
in the proton in spin, flavor, and momentum fraction, 
which have traditionally been ignored in analyses of dou- 
ble parton scattering, but only a small correlation with 
the transverse momentum kj^ . The uu and ud dPDFs are 
not simply related to each other, or to the single PDFs u 
and d, because of the spin-flavor correlations in the pro- 
ton quark model wavefunction in Eq. ([7]). The results in 
this paper provide quantitative results for these diparton 
correlations, which will help in the experimental analysis 
of double parton scattering at the LHC. 
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Appendix A 

We collect simplified expressions for the bag model 
wave function in momentum space and the functions <p n 
needed for the PY projection. Several of these results 
were already obtained in Ref. [23] . The Fourier trans- 
form of the wave function is 



* m (p)= /d^e^^x) 



2^fti? 3 / 2 



Sl(K)Xr, 



\/n 2 - sin 2 ft \S2(n)p-(TXr, 
where k = \p\R and 



(Al) 



si(k) 



1 



sin(«; — ft) sin(ft + ft) 



K — ft 



ft 



s 2 (k) = 2jo(ft)ji(K) - ^ s i( K ) ) 



(A2) 



and \m is defined in Eq. ([5]). For the unpolarized and 
longitudinally polarized single PDFs this leads to 



m 2 m ~ 2 (ft 2 



(s 2 + sj + 2 Sl s 2 p z 



I'm 7, 75 I'm = (-1) 



sin 2 ft) 

3/2 

2(ft 2 - sin 2 ft) 
x [s 2 + s 2 (l - 2p 2 L )+2s 1 s 2 j5 z ] . (A3) 
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|*J [GeV] \kj [GeV] 

FIG. 7. Comparison of the dPDF spin structures as functions of x\ or |kj_|, keeping the other variable fixed. The left panels 
show the uu dPDFs, and the right panels show the ud dPDFs. The SuSu, SuSu* , AuAd and SuSd 1 distributions are negative. 



For the transversely polarized PDF we need The functions </>„, used in the PY projection, are 



dp) = o / rsm — T n (v). 

1 K (ft 2 - sin 2 fi)» J v"- 1 n 



(A5) 



with 



m/ . / 1 — cos 2il \ /l sin2f2\ 

xT, * » ,T, j. A t M lT , T(«)=(n — « sm2 V - (- + — — cos2i, 

1 1 1 sin 2Q 1 - cos 2il , . . . 

= [sl+sl(l-p±)+2s 1 s 2 p z . A4) 2 2ft 2S2 2 



n 2 - sin 2 ft 
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